Universality of the Wigner time delay distribution 
for one-dimensional random potentials 
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We show that the distribution of the time delay for one-dimensional random potentials is universal 
in the high energy or weak disorder limit. Our analytical results are in excellent agreement with 
extensive numerical simulations carried out on samples whose sizes are large compared to the lo- 
calisation length (localised regime). The case of small samples is also discussed (ballistic regime). 
We provide a physical argument which explains in a quantitative way the origin of the exponential 
divergence of the moments. The occurence of a log-normal tail for finite size systems is analysed. 
Finally, we present exact results in the low energy limit which clearly show a departure from the 
universal behaviour. 
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The problem of quantum scattering by chaotic or disor- 
dered systems is encountered in many fields ranging from 
atomic or molecular physics as well as in the scattering 
of electromagnetic microwaves. Some properties of the 
scattering process are well captured through the concept 
of time delay. This quantity, which goes back to Eisenbud 
and Wigner , is related to the time spent in the inter- 
action region by a wave packet of energy peaked at E. It 
can be expressed in terms of the derivative of the 5* ma- 
trix with respect to the energy. In the context of chaotic 
scattering the approach based on random matrix theory 
(RMT) provides a statistical description of the time de- 
lays. This problem was first studied by a supersymmetric 
approach and in |l3| by using a statistical analysis. 
This latter work provides a derivation for the one chan- 
nel case for the different universality classes. Recently 
it served as a starting point for [|| where the N channel 
distribution is shown to be given by the Laguerre ensem- 
ble of RMT. In spite of its success, such a description by 
RMT is not entirely satisfactory, in particular it does not 
apply to strictly one-dimensional systems |l(| for which 
strong localisation effects occur. Furthermore it does not 
shed much light on the physical mechanisms which are re- 
sponsible for the universal distribution. In this work we 
explore another approach by considering the scattering 
by a one-dimensional random potential. In this case the 
existence of universal distributions was first conjectured 
in H on the basis of a comparative study of two different 
models. This was further supported by fl(f | where the 
random potential is still of a different kind. 

The purpose of this letter is to present a new derivation 
that accounts for the universality and also to provide a 
physical picture that explains the origin of the algebraic 
tail of the distribution in terms of resonances. Further 
details will be given elsewhere [ p2[ . To begin with, let 
us briefly recall the model. We consider the Schrodinger 
equation on the half line x ^ 0: 



— tp k (x) + V(x)ip k (x) = k 2 <ip k {x). (1) 



We assume that V(x) has its support on the interval [0, L] 
and impose the Dirichlet boundary condition ip k (0) = 0. 
Therefore, for x ^ L stationary scattering states of the 
form 



(2) 



represent the superposition of an incoming and a reflected 
plane wave. Since there is only backward scattering, 
the reflection coefficient e lS ^ is of unit modulus and 
the Wigner time delay takes the form r(k) = d ^ . 
Such a model with a random potential can be viewed as 
a model of a disordered sample connected to an infinite 
lead. Instead of using the invariant embedding method 
as in |l5|,[l4| or stochastic differential equations S, our 
starting point is to relate the time delay to the wave func- 
tion inside the sample. This may be achieved by using 
the identity 



d /' dip* dtp 
dx \ dx dE 



d 2 ^ \ 
dxdE J 



(3) 



By integration over [0, L] one gets the so-called Smith 
formula 



T{k) = ~j^ dx\i> k {x)\ 2 -^sm5{k). 



(4) 



It expresses the time delay as the sum of a dwell time 
H and a term that can be neglected in the high energy 
limit. Inside the sample, the wave function and its deriva- 
tive may be written in the form ip k (x) = Af sin0(x)e^ x ^ 
and ip/.(x) = kN cos 6 '(x)e^ x > . The normalisation factor 
|VV| = e~£( L ) is fixed by matching the wave function at 
x = L with the scattering states (d). We now consider 
the case where V(x) is a random potential. In this case 
the growth or decay of the envelope of the wave 

function is measured by the Lyapunov exponent 7 (in- 
verse localisation length A = 7 _1 ). In the high energy 
limit, the envelope is a slow variable, while the phase 



1 



6(x) presents rapid oscillations on a scale There- 
fore, in the high energy limit one can integrate over the 
fast variable in (0) and get: 



1 f 1 



da;e 2 «W-« L ». 



(5) 



This representation of the time delay holds for any real- 
ization of the disordered potential. It therefore captures 
all the statistical properties of r(fc) once the distribution 
of £(a;) is known. Denoting by x c the correlation length 
of V(x) and assuming that x c and are the small- 
est length scales of the system, then it was proven in |l]] 
that the variable is a Brownian motion of the form: 
£(x) = r yx + ■ } /jW(x) where W{x) is a normalised Wiener 
process ((W(x)) = 0, (W(x)W(x')) = min(x,x')). Thus 
the Lyapunov exponent 7 controls both the drift and the 
fluctuations. Using the scaling properties of the Brown- 
ian motion then gives the following identity in law: 



r(k) 



(law) 



due -2u+2W(u) _ 



(6) 



describe the evolution of the phase take a discrete form 
which can be implemented conveniently in a numerical 
simulation. The distribution of the time delay that we 
have obtained numerically in this way is in perfect agree- 
ment with equation (|^) as soon as the high energy regime 
is reached. For example we compare on figure [l] the an- 
alytical expression (Q) with the corresponding numerical 
result for a regime p -C v <C k. The simulation was car- 
ried out for a system with 10 5 impurities of weight v = 1, 
distributed with an averaged density p = 0.1. The energy 
considered corresponds to k = 10 which is related to a lo- 
calisation length A = ^ = 8000. The ratio L/X = 125 is 
sufficiently large for tfie limit distribution to be reached. 
The numerical calculation is based on a statistics of 50000 
values. It shows that the algebraic tail is well reproduced 
by (0) for 2000 times the typical value r typ = 200. Let 
us stress that there is no adjustable parameter to fit the 
numerics. The only parameter entering in the analytical 
expressions is the localisation length which is known for 
each kind of disorder. 



This representation of the time delay as an exponential 
functional of the Brownian motion p3| , ^8p |] allows to de- 
rive a number of interesting results: 
(?) existence of a limit distribution (r fixed, L — ► 00) with 
an algebraic tail |17) : 
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P(r) 



A 



2fcr 2 



(7) 
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(m) Linear divergence of the first moment and exponen- 
tial divergence of the higher moments [H : 



(r(fc)) 

(T(*r> 



L 

k 



(8) 
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FIG. 1. Time delay distribution in the localised regime 
L S> A. Comparison between the numerical calculation and 
expression (Q). In the inset: tail of the integrated distribution, 
numerical and analytical. 



In H we have shown that these results hold for two dif- 
ferent models of random potential in the localised regime 
(L>A). 

In order to test the analytical results in the above men- 
tioned regime it is convenient to choose a model suitable 
for numerical simulations. For this purpose we have con- 
sidered the case where the random potential is given by 
a sum of delta functions of the same weight v, randomly 
dropped on the half line with an average density p (the 
so-called Frish and Lloyd model)^. The equations that 



This model coincides with the Gaussian model (Halperin 
model) considered in |i| in the limit of a high density of im- 
purities (v <C k <C p). 



The derivation of the statistical properties of r given 
above allows one to understand the universality of the 
result but on the other hand does not shed much light on 
the physical mechanisms which are responsible for the oc- 
curence of an algebraic tail. In the following we propose 
a physical picture based on the existence of resonances 
that explains the leading exponential behaviour of the 
moments. The starting point is to realize that the reflec- 
tion of the incident wave on the random potential can 
in fact be viewed as a resonance tunneling process 
Indeed there exists a representation of the time delay as 
a superposition of resonances of energy E a and width T a 
in the form fll2|l 



2 



TJ2 



(E-E a ) 2 + T 2 JA 



(10) 



Obviously the dominant contribution r 



when E is in a window of width L n 



- is achieved 
centered at E a , and 
this will occur with probability ^ where A is the mean 
level spacing. In order to estimate the width we may 
assume that a discrete level E a localised at x$ will be 
broadened by its coupling to the continuum of states 
through the end point x = L. We may therefore set 
r ~ Q-2-f L {L-x ) -where ^ L } s ^ ne Lyapunov exponent in 
the finite system. Assuming that xq is uniformily dis- 
tributed on [0, L\ and decorrelated from 7/, one obtains 
the estimate: 



dzp 
L 



&1lp{il) 



r 1 



(11) 



Since £(£) = Ljl defined previously is a Gaussian pro- 
cess, the distribution of the finite size Lyapunov exponent 



v{il) 



One finally obtains 



2-777 



n\ p2n(n— 1)L/X 



(12) 



(13) 



A more refined derivation |2^] allows to recover the gross 
behaviour of the pre-exponential factor. This demon- 
strates that this particular behaviour of the moments has 
his origin both in the exponentially small widths of the 
resonances and in the fluctuations of the Lyapunov ex- 
ponent for the finite size sample. 

This physical picture shows that the problem of the 
time delay distribution related to the question of quan- 
tum relaxation in disordered systems. In this respect, the 
fact that P(t) is a broad distribution is rather puzzling. 
Indeed it was recently argued that at large time, the prob- 
ability distribution of various quantities, such as the con- 
ductance of one-dimensional disordered sample, exhibits 
a log-normal tail fl9|| . In order to clarify this point, in- 
stead of considering as before the regime r fixed L — > 00 
which leads to (|t]) , we have studied for fixed L the tail of 
the distribution in the limit t — * 00. In order to extract 
the asymptotic behaviour, it is convenient to consider the 
characteristic function <f>(p,L) = J a dr e~ 2kpT P(t; L) 
given in Jig] , If the conjugated variable p is chosen in 
a range ^e~ jL «p« 7e _v ^, the characteristic func- 
tion exhibits the following behaviour 



20Fe~TT 
(27L) 3 / 2 



In "f/p 



l + O 



e s ~i L . 



(14) 



which suggests the existence of a log-normal tail for the 
distribution 



P{t-L) ~exp-^ In 2 (2fc 7 r) 



(15) 



in the range e 7i ^> 2kjr ^> , Although we were 
not able to derive the behaviour of the distribution when 
e lL <C 2k-fT, the fact that the most divergent part of the 
moments grow like e 2n L l x suggests that the distribution 
is still log-normal. 

As an aside remark let us mention that a log-normal 
distribution of the time delay also occurs in the study of 
the random mass Dirac model at the critical point E = 
pl[ . There, the authors provide a representation of t 
which is similar to equation (^), except that the drift 
term is absent. This problem can also be analysed by 
using the approach given in 




FIG. 2. Second cumulant of the time delay in the ballistic 
regime L A. Comparison between numerical results and the 
analytical result (r 2 ) — (r) 2 = 6 ^ k S^_ pv y^ L 3 . the parameters 
are v = 0.001, p = 100 and k — 1. In the inset: time delay 
distribution for L = 100 (10 4 impurities). 



At this stage we have only considered a localised regime 
for which the size of the system is large compared to the 
localisation length L ^> A. Another interesting case is 
the ballistic one characterised by L <C A. In this situa- 
tion, the dimensionless variable 7L which arises in @ is 
small compared to 1 and the argument of the exponen- 
tial typically remains small compared to 1, which allows 
one to expand the exponential. The resulting expression 
for the time delay is given by a linear functional of a 
Gaussian quantity and has itself Gaussian fluctuations 
characterised by a first moment (r) = 4 and a second 
cumulant (t 2 ) — (r) 2 ~ ■^ I L i . We have checked numer- 
ically these results with the delta impurity model. We 
have considered a regime where it reproduces the high 
energy features of the Gaussian model: v <C k p 
and (k 2 — pv) ^s> (pv 2 ) 2 / 3 . In this regime one has to 
take into account the averaged value of the disorder pv 
and replace k by \/k 2 — pv in all previous expressions: 

the localisation length is thus given by A = B< " k p ~-i pv ^ 



3 



and the moments of r now read 



and 



(r 2 )-(rf 



pv 

6(k 2 — pv) 2 



\/k 2 -pv 

L 3 . In figure H we compare the nu 



merical result to the Gaussian distribution where the pa- 
rameters are given by the previous expressions. The cal- 
culation is performed for a ratio L/X ~ 1.4 10~ 3 . 10000 
values of r were calculated. 




FIG. 3. Time delay distribution in the localised regime 
L > A at low energy. 



Let us close the paper with some remarks on the low- 
energy regime. In this case one is more sensitive to the 
precise nature of the disorder, therefore one can guess 
that universality must break down. The distribution of 
the time delay will now depend on the nature of disor- 
der. As an illustration let us consider the delta impu- 
rity model; we predict |p2f , in the low density regime 
k -C p -C v, an exponential tail for the distribution: 



P(t) ~ fcpY(r - i) e- kp ^ T -^\ 
kv 



(16) 



where Y(x) denotes the Heaviside function. We have 
checked that this expression is in very good agreement 
with the numerical results. The numerical computation 
was performed for v = 1, p — 0.1, k = 0.01 for 1000 
impurities. The resulting distribution presented on fig- 
ure H is based on 50000 data sets. The behaviour at the 
origin is more subtle than the one given above, neverthe- 
less (|l6|) gives the correct scale on which the distribution 
vanishes at the origin. Another indication that universal- 
ity breaks down is the fact that within this same model, 
the low energy regime with a high density of impurities 
k -C v <C p leads to different distributions, though still 
characterised by an exponential tail. 
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